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A  new  ccnillrtueu^-fW.  direct-connect,  h]gh-enitialpy,supetsoiiR  eoiiibusticm  research  facility  is  described.  This 
lest  facility  provides  combustor  inlet  Sow  conditions come^pondmg  to  Bight  Mach  number*  between  3.5  and  7,  at 
dynamic  pressu  res  up  to  95.8  kPa.  Most  of  the  major  com  portents  of  the  new  facility  are  water  cooled  (including 
the  vitiated  healer,  the  instrumentation  and  transition  sections,  and  the  facility  norale  and  isolalorsh  the  current 
exception  is  the  variable-geometry  he  at- sink  combustor.  A  variety  of  conventional  and  advanced  instrumentation, 
including  a  steam  calorimeter  and  a  thrust  stand*  exists  for  accurate  documentation  of  combustor  inlet  and  exit 
conditions  and  performance  parameters.  In  a  recent  calibration  effort,  pitot  pressure  surveys*  total  temperature 
surveys,  and  wall  static  pressure  distributions  were  obtained  tor  a  wide  range  of  Inlet  conditions  LiSiny  Mach  1J* 
and  2J  facility  nozzles.  In  addition*  three-dimensional  numerical  simulations  of  each  test  case  were  completed. 
Results  from  the  computations  com  pare  favorably  with  experimental  results  for  ail  cases  and  yield  estimates  of  the 
integral  boundary-layer  properties  at  the  isolator  exit. 
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Introduction 

SI  MCE  1 995,  the  Ai  r  Force  Rcse  arch  Laboratory’s  Propul  sion  Di¬ 
rectorate,  Aerospace  Propulsion  Office  (AFRL/PRA)  has  been 
developing  a  modem,  well-characterized,  and  wdl-doeumcntcd 
direct-connect  test  facility  for  full-scale  scram  jet  combustor  de¬ 
velopment.  Facility  calibration  and  direct-connect  combustor  ex¬ 
periments  have  been  underway  in  this  facility  since  approximately 
mid- 1997. 1-7  The  facility  provides  the  Air  Force  with  the  oppor¬ 
tunity  10  develop  and  study  unique  concepts  in  supersonic  com¬ 
bustor  fuel  injection,  flame  holding,  ignition,  and  inlet-combustor 
isolation,  with  the  ultimate  goal  of  accelerating  the  development  of 
technologies  required  for  hydrocarbon-fueled  dual -mode  sc  ramjet 
propulsion  systems.  Potential  applications  for  these  systems  include 
rapid  theaier response,  global -reach aircraft,  and  affordable  access 
to  space. 

The  new  facility  was  motivated  by  several  factors.  Compared  with 
hydrogen -fueled  concepts,  the  database  for  hydrocarbon-fueled 
scramjetsis  relatively  limited.  Furthermore,  there  is  a  strong  need 
for  the  Air  Force  to  develop  in-house  facilities,  expertise,  and  qual¬ 
ified  personnel  in  the  areas  of  high-speed  propulsion  systems  test¬ 
ing  and  analysis.  The  goals  of"  the  current  effort  include  support of 
the  Hypersonic  Technology  {HyTech)  program"  (development  of  a 
Mach  4-S  hydrocarbon-fueledscramjet  missile)  and  exploration  of 
a  Uern  a  le  concepts  for  improve  ment  of  superson  iccomhu  stor  pert'or- 
mance.The  facility  is  designed  and  developed  to  perform  relatively 
long-duration*  steady-state direct-connect  combustor  tests.  Care  is 
being  taken  to  characterize  and  accurately  document  all  aspects  of 
this  facility.  Modem  measurement  techniques  and  Instrumentation 
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arc  being  incorporated  for  evaluation  of  scramjet  performance,  in¬ 
cluding  thrust-based,  chemical-based,  and  calorimclcr-bascd  com¬ 
bustion  efficiency. 

Facility  Descrip  Linn 

The  following  sections  describe  the  testing  and  simulation  capa¬ 
bilities  of  the  new  test  facility,  some  details  of  the  various  hardware 
components  that,  comprise  the  flowpath,  and  the  available  instru¬ 
mentation  and  control  systems.  The  test  cell  receives  continuous 
airflow  of  13.6  kg/s  at  5.17  MPa  and  a  maximum  temperature  of 
922  K  with  20.7  k  Pa  continuous  exhaust  from  the  Research  Air  Fa¬ 
cility.  Both  liquid-  and  gaseous-fuel  systems  are  a  vailablednciuding 
pumped  JF-4,  pressurized  JP-7.  ethylene,  and  hydrogen.  Uquid- 
and  gaseous-oxygen  systems  are  available  for  makeup  oxygen  in 
the  vitiated  heater.  A  recirculated  cooling- water  system  provides 
158  1/s  at  4-83  kPa;  raw  dump  water  at  2.41  MPa  is  also  avail¬ 
able.  The  entire  flowpath  is  secured  to  a  thrust  stand  for  direct  mea¬ 
surements  of  the  thrust  produced.  This  measurement  may  be  com¬ 
bined  with  wall  static  pressure  measurements  and  a  performance 
analysis  code  to  estimate  the  combustion  efficiency.  Additional ly„ 
the  energy  losses  through  the  various  water-cooled  components, 
coupled  wuth  temperature  measurements  from  a  steam  calorime¬ 
ter,  allow  calculation  of  combustion  efficiency?  Additional  details 
regarding  the  hardware  design  and  fabrication  may  be  obtained 
elsewhere.1** 

Testing  and  Simulation  Capabilities 

Johns  Hopkins  University/Applied  Physics  Laboratory's  (J  HLl/ 
APLFs)  Ramjet  Performance  Analysis  (RJPA)jUS  code's  in- 
let/diffuscr  model  was  used  to  relate  the  actual  flight  conditions  of 
a  hypersonic  vehicle  (and  the  corresponding  combustor  inlet  con¬ 
ditions)  to  tbe  test  cell  operational  capabilities  and  to  generate  the 
facility  map  shown  in  Fig.  1 .  This  map  indicates  the  range  oi‘  flight 
conditions  (altitude,  Mach  number  and  dynamic  pressure)  that  can 
be  simulated  in  the  test  cell  based  on  the  supply  air  pressure,  tem¬ 
perature,  and  mass  flow  constraints. 

The  inlet/diffuser  calculations  from  RJPA  w-cre  also  extended  to 
dc  sign  the  ex  perimentsconducted  In  the  testcell .  RJPA  was  executed 
in  an  iterati vemanner  such  that  the  follow  i  ng  two  cornel  ad  ons1 3  weir 
simultaneously  satisfied  for  the  range  of  flight  simulation  conditions 
presented  in  Fig.  1  (3.5  <  jV/.:,  <7  and  23.9  k Pa  <  Qti  <95.8  kPa). 
These  calculations  yielded  predictions  of  the  actual  flow  properties 


at  the  combustor  inlei  (T^mm  F4,aa>  and  AA.an )  at  each  flight 
condition  assuming  the  inlet  of  the  hypersonic  vehicle  varies  to 
satisfy  Eqs.  (1)  and  (2),  that  is,  die  inlet  is  optimized  for  each  flight 
Mach  number. 

(Ad/Aj)  =  -3.5  +  2.17*^0.017^  (1) 

(P4/P0)  -  —3.4  -  3.5*4  +  0.6_Wr;  (2) 

In  the  direct-connect  rest  facility,  four  facility  nozzles  (with  de¬ 
sign  Mach  numbers  of  3.8,  2F2,  2.6.  and  3.0)  are  available  to  es¬ 
tablish  the  conditions  entering  the  combustor.  If  the  vitiated  air¬ 
flow  through  the  nozzle  is  assumed  to  be  isentropic  with  a  specific 
heat  ratio  of  y4  =  1.3,  then  the  facility  stagnation  conditions  may  be 
computed  such  that  the  static  temperature  and  pressure  entering  the 
combustor  match  those  dcicrmined  using  RJPA  (Tt.^  =  and 
^4.sini  “  respectively).  Table  3  pre  sents  a  summary ofthe  sim¬ 
ulation  capabilities  using  the  various  facility  nozzles.  Tabic  3  also 
shows  how  well  the  combustor  Mach  number  and  velocity  obtained 
using  the  available  nozzles  compare  with  the  values  calculated  using 
RJPA.  The  condition  s  where  each  fad  iity  nozzle  ope  rate  son-de  sign  f 
(that  is,  actual  values  of  T4,  P. '4,  jW4i  and  fA  closely  match  the  RJPA 
output)  are  noted.  Because  each  facility  nozzle  can  exactly  match 
all  parameters at  only  one  condition* compromises  will  exisi  in  both 


Table  t  Facility  test  conditions  using  •y*  ■  1.3  and  forcing  =  r4[ Hd  and  /V.am  = 


JWfl 

£fl-kPa 
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A  U^% 
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Facility 


Fig,  3  Test  cell  22  hardware,  upstream. 


Mach  number  arid  velocity  if  the  nozzle  ls  operated  at  conditions 
away  from  that  design  point. 

Hard  wart 

Figure  2  is  a  schematic  of  the  facility  with  flow  from  left  to  right, 
Supply  air  enters  a  modified  M&rquardt  sudden  expansion  (SuE) 
vitiator  capable  of  sustaining  lemperat Lines  to  about  2500  K.  The 
vitiator  is  fitted  with  an  igniter  system  (Hj/air)  and  can  be  fueled 
with  both  liquid  and  gaseous  fuels.  For  the  tests  described  here, 
IP-4  was  used  to  fuel  she  vitiator.  A  water-cooled  instrumentation 
section  placed  downstream  of  the  vitiator  permits  stagnation  tem¬ 
perature  and/or  pressure  measurements  with  traversing  probes,  wall 
static  pressure^as  sampling  ports,  and  eight  thermocouple  ports  in 
the  inlet  and  outlet  flanges.  A  water-cooled  transition  flange  com¬ 
presses  the  vitiated  airflow  from  ax i sy mm etr ic ( 2-54 mm  i.d-lto  two 
dimcnsional(57.2  x  17  7,8  mm)  as  it  enters  the  facility  nozzle.  These 
components  appear  in  the  photograph  shown  in  Fig.  3. 

Four  interchangeable  water-cooled  facility  nozzles  are  available 
to  expand  the  air  to  supersonic  conditions,  simulating  desired  com¬ 
bustor!  nl  et  flow  cond  itions  .  The  nozzle  section  has  23  static  press  u re 
taps  on  its  side  wail  for  documenting  the  axial  static  pressure  pro¬ 
file  from  the  nozzle  inlet  to  exit.  Two  removable  isolator  sections 
downstream  of  the  nozzle  (water  cooled,  each  3G5  mm  long  with  a 
total  of  1 20  static  pressure  taps  on  top  and  bottom  walls)  are  used  to 
contain  the  pnecombustiort  pressure  rise.  The  isolatorsections  incor¬ 
porate  a  0r75-deg  divergence  along  ihc  bottom  wall  to  compensate 
for  boundary-layergrowth. 

Immediately  following  the  isolator  sections  is  a  variable- 
geometry  heat-sink  combustor  with  a  flexible  upper  wall  currently 
set  to  provide  an  area  rati  o  of  approximately  2.5.  The  comhustoral  so 


has  removable  in  sens  on  all  four  walls,  allowing  optical  access  and 
installation  of  a  variety  of  instrumentation.  The  replaceable  inserts 
provide  a  wide  parametriedcsign  space  for  fuel  injection  and  flame- 
holding  concepts.  The  initial  configuration  incorporates  four  Sow- 
angle  (15  de^)  flush- wail  fuel  injectors  upstream  of  a  cavity-based 
flameholder.-*  Liquid  and  gaseous  hydrocarbon  fuels  are  available 
for  the  combustor.  An  air  throttle  is  available  to  facilitate  combus¬ 
tor  ignition.  The  air  throttle  supplies  controlled  quantities  of  com¬ 
pressed  air  into  the  combustor  at  a  desired  axial  station.  Its  duration 
is  variable  {minimum  of  approximately!  s):  a  typical  cycle  involves 
between  4  and  10  s  of  airflow. 

A  calorimeter  instrumentation  section  that  houses  water  sprays, 
rakes,  and  probes  connects  the  combustor  to  a  calorimeter  instru¬ 
mented  with  thermocouple  sand  heat -flux  gauges.  The  calorimeter 
exit  connects  to  an  el  bow  through  which  the  flow  exits  the  test  cell  to 
an  exhaustersy  stem.  The  photograph  m  Fig.  4  shows  the  combustor, 
calorimeter  instrumentation  .and  calorimeter  sections. 

The  test  hardware  also  includes  a  spray-cooled  hcarbon-steel  cal¬ 
ibration  section  that  can  be  attached  to  the  exit  of  the  facility  nozzle 
or  either  isolator.  This  section  has  ports  for  water-cooled  travers¬ 
ing  probes,  allowing  span  wise  and  transverse  total  temperature  and 
pitot  pressure  surveys  at  the  nozzle  and/or  i  sola!  or  exit  planes.  This 
device  incorporates  a  sudden  area  change  at  the  probe  station  ^o 
that  the  probe  may  be  positioned  above  and  below  the  actual  top 
and  bottom  walls.  The  attachment  of  this  section  to  the  exit  of  either 
isolator  results  in  the  formation  of  a  weak  compression  wave  from 
the  bottom  wall  as  the  flow1  turns  through  0.75  deg  toward  the  rig 
cen  terl  me.  A  re  motel  y  acm  aied  butterfly  val  vc  in  stal  led  down  stream 
of  this  section  permits  studies  of  a  simulated  precombusiion  shock 
l  ruin  positioned  in  the  facility  isolator  sections. 
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Fig.  4  Test  cell  22  hardware,  downstream. 


Itistnmienrtartkm 

A  CAM  AC  process  1/Q  system  is  used  for  both  data  acquisition 
and  process  control.  l\vo  Sun  630MP  workstations  serve  as  the  data 
acquisition  and  facility  control  systems.  The  data  acquisition  system 
has  4  3  6  channels  of  analog  input,.  64  channels  of  digital  I/O,  and 
4f>  cha  n  ne  is  of  analog  output .  The  C  A  M  AC  crates  are  connected  to 
the  two  workstations  via  a  fiber  optic  small  computer  system  -SCSI 
interface.  A  Pressure  Systems,  Inc.  (PSI)  8400  pressure  scanning 
system  currently  consisting  of  400  channels  with  real-time  display 
and  data  reduction  is  also  available 

Total  mass  flow  rate  through  the  rig  at  each  test  condition  is 
estimated  using  two  independent  methods.  The  first  method  sums 
the  flow  rates  of  facility  air,  vitia  tor  fuel,  and  makeup  oxygen,  mea¬ 
sured  using  an  orifice  p  late  and  two  tu rbine  meters .  respective!  y.The 
second  method  uses  the  stagnation  pressure  and  temperature  mea¬ 
surements  in  the  instrumentation  section*  the  known  facility  noz¬ 
zle  throat  area,  and  simple  one-dimensional  gasdynamic  relations. 
The  close  correspondence- between  these  two  estimates,  generally 
within  1-39**  inspires  confidence  in  the  quality  and  reliability  of 
i  ns  trumen  tation  and  an  al  ysis  i  n  th  e  test  cell  3  tags)  ation.  temperatu  re 
is  currently  measured  using  thermocouples  (for  temperatures  up  to 
around  1 200  K). 

Remote  monitoringof  test  activities  is  providedby  five  videocam- 
era.H  placed  at  strategic  locations  throughout  the  test  cell.  Four  cam¬ 
eras  placed  around  the  thrust  stand  monitored  the  facility  hardware. 
A  hand- held  Sony  ft-mm  camera  monitored  the  traversing  probe. 

Calibration  Efforts 

Experimental  and  numerical  studies-  were  conducted  to  calibrate 
two  facility  nozzles  (Mach  l.-R  and  2.2)  in  conjunction  with  the 
isolator  sections.  Table  2  presents  the  stagnation  conditions  corre¬ 
sponding  to  the  12  calibration  cases  examined.  The  following  sec¬ 
tions  describe  the  expcrimcnla  I  and  computational  methods  used  to 
diagnose  the  flow-field  and  the  results  of  the  cal  ibral  ion  studies.  For 
l>rev  try,  resu  I  ts  from  on  e.  case  for  each  fae  i  I  ity  n  ozzlc-  (c  a  ses  3  and  9 ) 
are  presented  m  this  paper. 

Experimental  Details 

Two  separate  water-cooled  traversing  probes  were  used  for  hard¬ 
ware  calibration.  The  first  probe  contained  eight  pitot  pressure  tips 
covering  approximately  25  mm  in  the  span  wise,  direction.  The  sec¬ 
ond  probe  consisted  of  four  p  i  tot  pressure  t  ips  an  d  two  total  tempc  ra- 
t  ure  tips  over  the  same  span .  The  pilot  pres  sure  tip  s  were  made  from 
hypodermic  tubing  (0.53  mm  i,d.).  Total  temperatures  were  mea- 


TabkL  2  Experimental  and  cum  pula  tiuna  L 
calibration  test  conditions 


Gas* 

^U.-wCI 

Ph  kPa 

1 

1.8 

833 

372 

0.97, 2  78. 3.52 

2 

18 

778 

703 

0.94,2.76,3.54 

3 

1.8 

972 

372 

0.96,2.72,3.51 

4 

1.8 

917 

752 

0.94,2.84,3.49 

5 

I.B 

1139 

365 

0.96. 2.59,  3.34 

6 

LB 

1083 

745 

0.94, 2.51,3-50 

7 

2.2 

902 

689 

0,94. 2.55,4  46 

8 

2.2 

917 

1351 

0.92,3.12,4.05 

9 

2.2 

1083 

710 

0  94. 3.00, 4  50 

10 

2.2 

1063 

1441 

0.9:2,3.05, 3.86 

11 

2-2 

1259 

703 

0.93*  2.65*4.35 

12 

2.2 

1232 

1420 

0.92,3.64,  3.93 

sured  using  fine  gauge  thermocouples  set  inside  aspirated  ceramic 
diffusers.14  Both  probes  were  inserted  into  the  calibration  section 
from  the  top  and  were  connected  to  a  computer-  controlledelecLronic 
actuator  for  remote  positioning.  Press  ures  were  measured  using  the 
PSI  system  while  the  thermocouple  signals  were  recorded  via  the 
CAM  AC  data  acquisition  system. 

Co  m  pulationa  [  IX  a  i  Is 

Three-dimensional  simulations  were  performed  at  low  backpres¬ 
sure  for  each  of  the  1 2  eases  outlined  in  Table  2.  These  computations 
provide  resul  ts  for  com  pari  son  w  ith  the  experimen  ta!  measu  re  men  t  s . 
The  computational  results  provide  a  suitable  means  for  calculating 
t  he  bound  ary-  lave  r  properties  (bound  a  ry-lay e  r  th  ickncss  & ,  the  d  is  ■ 
placement  thickness 


and  the  momentum  thickness 


at  the  isolator  exit  plane.  Additional  simulations  were  performed 
at  elevated  backpressures  for  cases  1  and  9.  The  method  used  to 
calculate  the  flowficlds  at  low  back  pressure  will  first  be  described 
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followed  by  modifications  employed  for  the  elevated  back-pressure 
simulations. 

The  computational  domain  included  the  subsonic  and  supersonic 
portions  of  the  nozzle,  the  isolator,  and  die  47. 6-inm-long  constant 
area  calibration  section  that  was  attached  at  the  end  of  the  isolator. 
Span  wise  symmetry  was  assumed;  hence,  the  domain  extended 
from  the  tunnel  centerline  to  the  side  wall  in  the  spanwise  direc¬ 
tion.  The  calculations  were  performed  with  the  VULCAN  Navier- 
S  tokes  code , 1 5  a  derivative  of  the  LARCK  code  devel  oped  at  NASA 
Langley  Research  Center-  The  calculations  employed  the  low- 
diffusion  f\u%  vector  split  scheme  of  Edwards'6  for  the  evaluation 
of  the  inviscid  fiuses.  The  MUSCL  parameter  was  set  to  one-third 
to  minimize  truncation  error.  The  i.hin-3ayer  viscous  terms  were 
computed  in  the  three  coordinate  directions.  ThjtmSence  was  mod¬ 
eled  with  the  Menlcr17  two-equation  BSL  turbulence  model.  At 
solid  surfaces,  the  wall-matching  procedure  of  Wilcox15  was  em¬ 
ployed  to  reduce  grid  requirements.  The  turbulent  Prandtl  num¬ 
ber  was  set  to  0.89.  The  subsonic  portion  of  the  nozzle  extended 
a  short  distance  downstream  of  the  throat  and  was  solved  cllipti- 
tally,  whereas  the  sections  downstream  were  solved  with  the  space 
marching  algorithm, 

The  simulated  vitiated  a irslream  composition  was  obtained  from 
the  experimental  flow  rate  measurements  assuming  complete  com¬ 
bustion.  The  resulting  mass  fraction s of  N2,  O-,,  CG2,  and  H+Q  were 
0.722,  0.232,  0.033,  and  0  0 13  for  case  1  and  0.688, 0.2 32,  0.05  S> 
and  0.022  for  case  0,  respectively.  The  calculations  were  performed 
with  a  single  species  that  had  the  molecular  weight  and  thermody¬ 
namic  coefficients  set  equal  lo  the  vitiated  heated  air  mixture  valucs- 

The  inflow  boundary  was  prescribcdby  fixing  the  stagnation  pres¬ 
sure  and  temperature.  The  stagnation  temperature  at  the  inflow  of 
the  nozzle  was  obtained  by  adjusting  the  measured  temperature  at 
the  exit  of  the  vitiated  heater  for  the  heat  loss  through  the  instru¬ 
mentation  and  transition  sections.  The  nozzle  inflow  total  temper¬ 
atures  for  cases  I  and  9  were  SIS  and  1058  K,  respectively.  The 
outflow  boundary1  condition  for  the  elliptic  region  was  first-order 
extrapolation.  The  top,  bottom,  and  side  walls  employed  a  no-slip 
iso  thermal  boundary  condition.  The  wall  temperature  in  the  noz¬ 
zle  and  isolator  sections  was  determined  by  iteratively  performing 
a  series  of  two-dimensional  calculations  until  the  calculated  heat 
loss  matched  the  cxpcrimcnlally  measured  heat  loss.  These  wall 
temperatures  were  then  imposed  on  the  three- dimensional  calcula¬ 
tions.  The  measured  heat  loss  exceeded  the  heat  loss  resulting  from 
each  three-dimensional  calculation  by  Less  titan  5%.  The.  calibra¬ 
tion  section  was  cooled  to  an  unknown  extent.  Ln  these  simulations, 
the  calculated  heat  Loss  through  this  section  was  slightly  overpre- 
dictccL  thereby  compensating,  to  some  extent,  the  underpredicted 
heat  loss  in  I  he  nozzle  and  isolator  sections,  The  resulting  nozzle 
and  isolator  wall  temperatures  were  400  and  330  K.  respectively, 
for  case  1,  and  455  and  380  K,  respectively,  for  case  9  Finally,  a 
symmetry  boundary  condition  (zero  span  wise  velocity  component 
and  extrapolation  ot'  all  other  flow  quantities)  was  employed  along 
the  centerline  plane. 

The  calculations  were  advanced  in  time  with  a  diagonalized 
approximate  factorization  scheme19  with  the  Cou ran t-Frie dr iclis- 
Lewy  (CFL )  number  s  c  t  cqu  al  to  4 .0  and  5 . 0  for  the  elliptic  and  space 
marc  hi  n  g  region  s,  respecti  vely. The  con  v  ergencccritcrion  employed 
was  a  relative  error  of  the  L  2  norm  of  the  residua]  of  10-6.  The  mass 
flow  conservation  was  excellent  l  mass  loss  did  not  exceed  Q.3*&. 

A  sensitivity  study  was  conducted  for  case  l  using  three  differ¬ 
ent  grids.  The  coarsest  grid  had  171  points  in  the  axial  direction, 
40  points  in  the  Iran  s  verse  direction,  and  22  points  in  the  span  wise 
direction  (173  x  40  *  22).  The  finest  grid  doubled  the  number  of 
points  in  each  direction,  that  is,  341  *  79  *  43,  whereas  the  inter¬ 
mediate  gri  don  ly  doubled  t  he  n umberof  poi  ntsin  the  ax  ialan  d  tra  n  s- 
verse  directions, that  is,  341  *79  h  22.  Increasing  the  grid  density 
in  the  ax  ial  diree  lion  p  ro  vided  bet tc  rrc  solu  t  Ion  of  the  weak  p  re  ssu  re 
waves  in  the  nozzle  and  isolators,  whereas  coarse  resolution  In  the 
transverse  direction  yielded  thicker  boundary  layers-  The  solutions 
obtained  using  the  intermediate  and  finest  grids  were  examined  on 
the  rig  centerline.  Results  were  found  to  be  nearly  indistinguishable- 
Consequently,  the  calculations  for  the  remaining  cases  employed 
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Fig,  5  Computational  grid  for  Mach  1-8  fad  lily  nozzle  (ever)-  other 
point  removed]. 


Pic.  6  Mach  number  contours  from  Mach  1.8  fucilalv  nozzle  o  hala¬ 
tion. 


the  intermediate  grid  to  reduce  computational  time.  Figure  5  shows 
the  computational  grid  along  the  rig  centerline  for  the  Mach  1.8 
facility  nozzle  (every  other  point  has  been  removed  for  clarity).  The 
calculated  Mach  number  contours  along  the  rig  centerline  for  case  1 
appear  in  Fig.  6.  This  plot  indicates  that  the  nozzle  performs  as 
designed  with  an  exit  Mach  number  near  1 .8. 

Additionally,  calculations  were  performed  with  elevated  back¬ 
pressures  at  the  exit  pl  ane  of  the  calibration  section  for  ca  ses  I  and  9, 
The  isolator  and  calibration  sections  were  solved  together  eJlipti- 
cally  with  the  flow  variables  at  the  inflow  boundary  set  equal  to  the 
calculated  va|u  es  at  the exitpl anc  of  the  nozzle ,  and  a  11  flow  v  ariablcs 
at  tb  e  ou  tflo  w  bou ndary extiapolatcdfrom  the  i  ntcrio  r  (except  for  the 
pressure,  which  was  imposed).  Only  results  from  the  highest  back¬ 
pressure  condition  will  be  presented  because  a  significant  region  of 
reverse  flow  was  present  at  the  outflow  plane  for  the  intermediate 
backpressure  condition.  At  the  highest  backpressure  condition  for 
case  9,  reverse  flow  was  present  in  a  small  region  in  the  comer  of  the 
duct  that  represented  approximately  0.03  and  0.40  &  of  the  exit  area 
on  the  coarse  and  fine  grids,  res  pc  ctivcly.Th  is  shou  Id  h  ave  a  min  ima l 
impact  on  the  results,  The  pressure  rise  for  the  highest  backpressure 
condition  for  case  1  extended  into  the  nozzle,  and  consequently, 
the  section  of  the  nozzle  downstream  of  the  throat  was  included  in 
the  elliptic  solution  domain.  The  wall  temperature  was  assumed  to 
be  the  same  as  in  the  non backp re ssured simulations. The  thin -layer 
approximation  for  the  viscous  terms  w!as  not  employed  in  the  sim¬ 
ulations  with  an  imposed  backpressure.  Solutions  were  obtained  on 
the  fine  grid  (34  3x  7 9  *  43)  that  had  additional  points  in  the  span- 
wise  direction  to  resolve  the  recirculation  zone  along  the  side  wall. 
Solutions  were  also  obtained  on  a  coarse  grid  with  every  other  grid 
pointremovedineachcoordinatedireclEonCl  TI  x  40  *  22)  to  inves¬ 
tigate  the  grid  sensitivity  of  the  results.  The  L2  norm  of  the  residual 
in  the  solutions  with  an  Imposed  backpressure  dropped  approxi¬ 
mately  three  orders  of  magnitude.  Convergence,  was  assessed,  prac¬ 
tically ,  by  monitoring  changes  in  mass- weigh  ted  one-dimensional 
quantities.  Convergence  was  considered  adequate  when  mass- 
weighted  one-dimensional  quantities  did  not  change  after  5000 
iterations. 

Mach  l.S  Facility  Nozzle  Results 

Figures  7- 11  present  the  results  ofthe  Mach  1  8  nozzle  calibra¬ 
tions  for  case  1 .  Walt  static  pressures  through  the  nozzle  and  isolator 
sections  appear  in  Fig.  7.  Included  in  this  plol  are  all  of  the  mea¬ 
sured  pressures,  the  data  obtained  from  the  numerical  simulation, 
and  a  theoretical  pressure  distribution  based  on  a  one- dimensional 
isentropic  analysis.  The  axial  position  is  sealed  by  the  nozzle  exit 
height,  and  x/H  =  0  corresponds  to  the  facility  nozzle  exit  plane. 
The  pressure  data  are  normalized  by  the  nozzle  exit  pressure  Pfl. 
With  low  backpressure,  static  pressure  profiles  through  the  nozzle 
and  isolator  sections  m  a  Lch  the  is  en  tropic  prediction  and  the  numeri¬ 
cal  simulation.  As  the  downstream  valve  is  closed,  the  pressure  at  the 
isolator  exit  rises,  resulting  in  a  system  of  shock  waves.  This  shock 
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Fig.  2  Normalized  wall  static  pressure  distributions  with  Mach  1-3 
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Fig.  9  Calculated  normalized  stalk  pressure  contours  from  case  1  at 
Fa JPA  =  3 . 5  2  i  n  the  vicinity  of  the  shock  Ira  in . 


Fig-  II  Normalized  isolator  exit  piano  total  temperature  distributions 
from  case  1  at  P^JPa  ~  0,97. 
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system  decelerates  the  approaching  flow,  Figure  7a  contains  mea¬ 
sured  static  pressure  distributions  for  two  levels  of  increased  back¬ 
pressure,  P.i/P4=2.7a  and  3-52.  The  pressures  along  the  top  and 
bottom  walls  of  the  isolator  sections  are  very  similar  in  each  profile. 
As  expected  .increasedexit  pressu  re  forces  the  shook  train  upstreani- 
Thes e  distribu  ti  onsillustratethc  effec  t  iven  cs  sof  the  isol  ators  in  con  - 
taining  a  combustorpres5ureri.se  thatiii  nearly  equivalentto  a  normal 
shock. 

The  computational  fluid  dyn  ami.es  (CFD)  prediction  at  the  higher 
backpressure  condition  are  shown  in  Fig.  7b  along  with  the  mea¬ 
sured  pressure  data  from  this  condition.  Along  the  centerline  plane, 
the  pres  su  re  ri  from  the  shock  train  begins  at  x  /  H  =  0,  which  i  s  ap¬ 
pro*!  mately  1 .5  duct  heights  upstream  of  the  measured  location  for 
the  beginning  of  the  pressure  rise.  The  flow  separates  initially  in  the 
comer  between  the  side  wall  and  the  bottom  walk  the  corresponding 
oblique  shock  causes  the  pressure  rise  to  occur  farther  upstream  on 
the  side  wall  than  along  the  centerline  plane.  Downstream  of  the 
nearly  normal  shock  along  the  centerline  plane,  there  is  a  series  of 
expansion  and  compression  waves  causing  the  pressure  to  decrease 
and  increase.  A  steady  solution  was  obtained  computationally.  The 
monotonic  rise  in  the  measured  wall  pressures  downstream  of  the 
shock  maybe  due  to  flow  unsteadiness  and  a  consequent  averaging 
of  the  crests  and  troughs  in  the  pressure  values,  The  shock  position 
and  overall  character  of  the  flow  is  represented  well  by  the  coarse 
grid  solution,  although  the  flow  quantities  are  more  diffuse  due  to 
numeric  a3  dissipation  - 

Figurc  5  shows  the  Mach  number  contours  along  the  centerline 
plane  and  several  crossflow  planes  for  the  solution  obtained  on  the 
fine  grid.  The  contour  of  zero  stream  wise  velocity  is  overlaid  in 
white  to  identify  the  boundary  of  the  recirculation  zone,  which  em¬ 
anates  from  the  comer  of  the  side  wall  and  bottom  wall  and  ex¬ 
tends  11-0 duct  heights  from  -1.2  <x/H  <9.8.  For  comparison, 
the  recirculation  zone  in  the  coarse  grid  solution  extended  9.4  duct 
heights  (—  G.9  < x/H  <  5.5).  At  the  exit  of  the  domain  the  Mach 
number  distribution  retains  a  significant  spanwise  variation  with 
lower  momentum  flow  toward  the  side  wall.  The  domain  extending 
from  - 1 .4  <  x/H  <3.7  is  outlined  in  black.  The  pressure  contours 
for  this  region  appear  in  Fig  9.  The  peak  pressure  ratio  along  the 
centerline  plane  is  P#/ ft  ^  3,57,  which  corresponds  to  97“%  ofthe 
normal  shock  pressure  rise  as  determined  from  the  local  values  of 
the  pressure.  Mach  number,  and  ratio  of  specific  heats  for  the  flow' 
upstream  of  the  shock.  The  footprint  of  the  oblique  shock  gener¬ 
ated  by  the  flow  disturbance  from  the  recirculation  zone  formed  in 
the  duct  comer  is  observed  on  the  bottom  wall.  The  pressure  rise 
along  the  bottom  and  side  wall  extends  slightly  over  one  duct  height 
upstream  of  the  shock  position  along  the  tunnel  centerline,  The  re¬ 
flection  of  the  oblique  shock  off  the  top  wall  of  the  duct  can  also  be 
seen  at  the  exit  plane  of  the  nozzle.  While  the  flow  in  this  region  is 
highly  three  dimensional,  the  pressure  held  downstream  of  this  re¬ 
gion  is  essentially  one  dimensional,  as  observed  io  the  downstream 
crossflow  plane  in  Fig.  9. 

Figure  10  shows  the  pilot  probe  mcasurcmcnls  obtained  just 
down  s  tream  of  the  aft  i  solatorcx  it  for  each  of  the  three  ba  ckpiessuie 
]  evels  exami  n  ed.  In  Fi  gs.  I  Oa  I  Oe  the  Iran  s verse  pos  irion  i  s  scaled 
by  the  facility  nozzle  exit  height,  and  y/H  =  0  corresponds  to  the 
top  wall  ofthe  isolator.  The  pi  tot  pressures  are  normalized  hy  the 
facility  stagnation  pressure.  Figure  10a  presents  the  experimental 
and  numerical  results  obtained  from  case  1  at  low  backpressure. 
The  eight  measured  pressures  are  very  -consistent  across  the  duct 
height  for  alt  cases  examined.  The  computed  profile  compares  well 
with  the  measured  values.  Figures  1  Ob  and  1  Oc  contain  normalized 
pitot  pressures  at  fti/ft  =  2.75  and  3.52,  respectively.  These  pro¬ 
files  again  exhibit  a  high  degree  of  consistency  from  one  probe  Lip 
to  another.  They  also  illustrate  the  effect  of  the  shock  train  on  the 
overall  flow  field.  At  the  highest  backpressure  level,  the  flow  exit¬ 
ing  the  aft  isolator  is  purely  subsonic.  The  CFD  results  shown  in 
Fig.  10c  suggest  that  the  shock  train  as  somewhat  asymmetric  from 
top  to  bottom.  Shock  train  unsteadiness  is  one  possible  cause  for 
these  differences.  The  experimental  data  represent  averages  over  a 
relatively  long  time  (and  potentially  many  oscillation  cycles),  and 
the  CFD  solution  is  steady  in  nature. 


Total  temperature  measurements  appear  in  Fig,  11  for  case  I 
operated  at  the  low  back  pressurecondition.  CFD  results  arc  overlaid 
on  the  measurements  and  compare  very  well  across  the  duct  height. 
This  plot  suggests  that  the  core  flow  is  not  strongly  affected  by 
Ihe  heat  loss  through  the  water-cooled  sections  (Tvt  for  case  1  is 
833  K). 

Mach  1,2  Facility  Nozzle  Results 

Figures  12-14  contain  results  of  the  Mach  2.2  nozzle  calibrations 
for  case  9  (refer  to  Tabic  2),  The  measured*  calculated,  and  isen- 
tropic  wall  static  pressure  distributions  shown  in  Fig.  1 2a  compare 
very  favorably.  As  with  the  Mach  3.8  results  (Fig.  7),  increased 


x/H 


■S  -4  0  4  B  12 


x/H 

b)  Comparison  of  CFI>  coarse  and  hue  grid  solutions  for  case  1 
tPJPA*4M) 


FLg.,  12  Normalized  wait  static  pressure  distributions  with  Mach  2.2 
facility  nozzle. 


Rig.  13  Normalized  isolator  ml  plane  pitot  pressure  distributions 
from  case  9  at  ft  \ /ft  »  &.  94 . 
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Table  3  Inl^ra  L  buunda  ry -layer  properl  Jest  a  I  a  ft  i  mj  i  a  Lcit  pa  J  t  pbme 


Case- 

Unit  Reynold? 
nu .,  m~ 1 

mm 

4-*S- 

mm 

&T?  6^;, 

mm. 

Sr/ij, 

h/r. 

*r/*r. 

&bJ$b 

rT  /Sr, 

1 

12-1  x  10s 

12.7, 12.3 

1.62,1,60 

1,35,1.29 

7.6, 7,7 

94, 9.5 

1.2, 1.2 

2 

25.3  x  L0* 

12  0, 11.4 

3.57,1.50 

149,1.12 

7.6, 7.6 

10.1,1 0.2 

1.3, 1.3 

3 

10.2  x  10s 

13.1,  12.6 

1.52,1.52 

1.42, 1.37 

8.6,  6.3 

9.2, 9.2 

1.1, 1.1 

4 

22.0x  U? 

12.3, 1LS 

!  .47, 1 .42 

1-27  J. 19 

8.4,  8.3 

9.7, 9-9 

1.2, 1,2 

5 

8.2  x  10s 

13.4,  13.1 

1.40,1  37 

1  55.1.50 

9.6,9  6 

8.6,3  7 

0.9, 0.9 

6 

17.2  x  10s 

12.5,  12  1 

1.37,1.35 

1. 37,1.29 

9.1, 9.0 

91,9.4 

1.0, 1.0 

7 

] (S_4  y  Iff5 

12.2,  12.2 

1.35,1.93 

1.22, 1.19 

6.6, 6.3 

10.0, 10.3 

1-5,1  62 

B 

3L2*  10* 

11-6,  11.5 

1.73,1.78 

L 12. 1.09 

6.7, 6-5 

104,10.6 

1.5,  l  6 

9 

13.1  x  10* 

1 2.6  J  2.6 

1  .S3.  L  88 

1,32,1.30 

69,67 

9.5,9  7 

1.4,  t  4 

to 

26.9  x  10* 

11.9  >11.9 

1.68, 1.75 

1.17,1.84 

7. 1 , 6.3 

10.2,10.4 

1,4, 1.5 

11 

10.8  x  10* 

1 3.2. 13.3 

1.68. 1.75 

1.42, 1.42 

7.8, 7.6 

9-3, 9.4 

1.2,  8 .2 

E2 

22.3  x  |0* 

1 2-3  J  2.3 

1.57,  ].62 

i  .27,1.24 

7,8,  7.6 

9-7, 9.9 

1,2,  1,3 

1200 

105'0 

900 


m 


450 
200 

0.0  0.2  0.4  0.6  0.S  3.0  1.2 

y/H 

Fiy.  14  Normalised  isolator  exit  plane  lota  1  tem  perature  distributions 
from  case  9  at  P^P*  ■  0.94. 

backpressure  results  in  the  formation  of  a  system  of  shock  waves 
positioned  in  the  isolator  sections.  Increasing  backpressure  forces 
this  disturbance  upstream,  At  this  in]c3  Mach  number,  the  isolator 
sections  effectively  contain  a  pres-sure  ratio  of  Pai/P^=4.50. 

The  CFD  pressure  data  at  the  higher  backpressure  condition  are 
shown  in  Fig.  L  2b  along  with  the  measured  pressure  data  from  this 
condition- There  is  good  agreement  between  the  CFD  and  measured 
data  regarding  the  position  of  the  pressure  rise,,  although  a  small  grid 
sensitivity  is  observed.  As  with  the  Mach  1  .3  nozzle,  the  pressure 
distribution  in  the  coarse  grid  solution  is  more  diffuse  than  in  the. 
line  grid  solution,  whereas  the  measurements  exhibit  a  monotonic 
pressure  rise. 

Figure  E  3  shows  the  measured  and  computed  pitot  pressure  dis¬ 
tributions  from  just  downstream  of  the  aft  isolator  exit  for  case  9 
operated  at  low  backpressure.  For  these  measurements,  the  probe 
with  four  pitot  tubes  was  used.  All  four  sets  of  measurements  agree 
closely  across  the  entire  duct  height.  Also,  the  numerical  results 
reproduce  the  features  of  the  measured  profile  very  well,  including 
the  weak  compression  wave  that  elevates  the  pitol  pressure  near  the 
bottom  wall  (0.7  <  y/fi  <0r9).  The  pitot  pressure  distributions  in 
ihe  lop  and  bottom  wall  boundary  layers  arc  well  represented  by  the 
CFD  results. 

Finally,  Fig.  14  shows  the  measured  and  calculated  total  tempera¬ 
ture  distributionsfrom  case  9  ope  rated  at  Low  backpressure.  As  in  the 
Mach  1 .8  nozzle  results,  the  measured  and  computed  total  temper¬ 
ature  distributions  compare  very  favorably  across  the  duel  height 
Again,  the  water-cooled  hardware  does  not  strongly  influence  the 
total  temperature  of  th=  core  flow  (7^  for  case  9  is  108-3  K). 

]n  WgraJ  Belinda  ry -Layer  F  ropertles 

Based  on  ihe  good  agreement  between  the  CFD  results  and  mea¬ 
sured  quantities,  the  calculated  profiles  were  used  to  determine  the 
integral  boundary-layer  properties  at  the  isolator  exit  plane  along 
l he  rig  centerline  for  each  of  the  1 2  cases.  The  trapezoidal  method 


was  use d  for  nume  rical  integ  rati  on .  Ve  locity  and  den  sity  val  ues  were 
obtained  at  the  cell  center  location*  and  the  grid  coord  inales  were  ob¬ 
tained  from  the  computational  grid.  The  boundary- layer  properties 
were  found  to  be  quite  sensitive  to  the  choice  of  the  edge  velocity 
Ut.  To  minimize  any  subjectivity,  the  average  velocity  in  the  region 
of  the  core  How  where  dU/dy  <  500  s”1  was.  defined  as  and 
the  edge  velocity  was  defined  as  Ue  =  0.99  x  U^,.  The  re  suits  of  the 
integration  process  arc  tabulated  in  Table  3. 

The  trends  m  the  calculated  values  for  5*  and  9  can  largely  be 
ex  pla  Ined  by  cons  i  deri  ng  the  di  srri  bution of  pf  pt  ■  For  a  given  di  s  trl- 
buti  on  of  U /  Ue ,  an  i  nc  re  ase  in  p  fp£  reduces*?* ,  the  effective  di  starsc  e 
occupied  by  ihe  boundary  layer,  and  Increases  6,  the  associated  drag 
force  on  a  flat  plate.  In  the  simulations  using  a  fixed  Mach  number 
nozzle,  p-fpt  increases  as  the  total  temperature  increases  due  to  the 
Corresponding  increase  in  heat  transfer  to  the  walls  of  the  water- 
cooled  hardware.  Con sequent  ly^  larger  values  are  obtained  for  6  as 
the  total  te mperatu  re  i  notes. ses .  The  uni  l  R  cy  nol  ds  nu  mbcrdccrca  ses 
with  increasing  total  temperature,  resulting  in  an  increase  in  &  and 
a  decrease  in  U j Ut  at  a  fixed  distance  from  the  wall  Finally,  the 
Increase  in  p/pe  is  larger  than  the  reduction  in  £//£/,.  resulting  in  a 
reduction  in  &  m  with  increased  total  temperature. 

Summary 

A  newly  developed  supersonic  combustion  research  facility  has 
been  described.  This  test  facility  offers  a  broad  flight  Mach  number 
and  dynamic  pressure  simulation  space  using  four  interchangeable 
facility  nozzles.  The  major  components  of  the  rig  are  water-cooled 
for  exposure  lo  relatively  long  test  limes  to  allow  for  steady -state 
measurements  of  combustor  performance.  Modem  data  acquisition 
and  control  systems  are  available,  along  with  a  variety  of  optical  di¬ 
agnostic  techniques  for  documentation  of  crucial  combustor  perfor¬ 
mance  paramctc  rs  i  itcludmg  the  i  nlet  co  n  di  Lion  s,  th  rust,  combustion 
efficiency,  combustor  he  at  loss,  and  wall  static  pressures.  The  results 
of  a  recent  effort  to  calibrate  two  of  the  facility  nozzles  demonstrate 
that  the  nozzle  and  isolator  sections  perform  as  expected.  Wall  static 
pressure  distributions  show  the  effectiveness  of  the  i  solatorsetti  ons 
at  conlaining  a  simulated  precombustion  shock  train.  CFD  results 
compare  well  with  the  wall  and  in-stream  measurements.  These  re¬ 
sults  have  been  a  na  ly  zed  to  obta  in  integral  bo undary-1  ayerpropertie  s 
at  the  aft  isolator  exit  plane. 
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